Phase space localization of antisymmetric functions 
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Upper and lower bounds are written down for the minimum information entropy in phase space 
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the wave function is restricted to belong to any of the proper subspaces of the Fourier transform 
operator. 
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1. Heisenberg uncertainty relations between position 
and momentum can be viewed as a constraint on the 
maximum allowable localization of quantum states in 
phase space, thereby displaying a bàsic principle of Quan- 
tum Mechanics. In this approach the localization is 
bascd on the sizc of the state. An alternative approach 
§, 1 1, |, | | @, § g [ÏÏJ 0, g H is based on using 
the information entropy of the state as a measure of its 
localization. To be specific, if ip( x ) is a normalized wave 
function, the information entropy of the corresponding 
state in position space is defined as 



s x W 



d d x\ip{ X )\ 2 \o g mx)\ 2 ), 



(1) 



(where d denotes the dimension of the space). The 
more localized the state the lower the entropy. For 
instance, for a one-dimcnsional Gaussian ip( x ) — 
(7T0 2 ) -1 / 4 exp(— x 2 /2a 2 ) (with a > 0), the behavior of the 
entropy S x = (1 + log(™ 2 ))/2 is qualitatively similar to 
that of the dispersion Ax — aj\[2. However, for a state 
composed of two well separated Gaussian functions of 
given widths, the dispersion would increase with the sep- 
aration whereas the entropy would not. This illustrates 
the different concept of localization in both approaches. 
The entropy in momentum space is defined similarly 

as 
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where 



$(fc) = / d d x 



ikx/h 



(3) 



For the one-dimcnsional Gaussian this gives ip(k) = 
(47ra 2 ) 1 / 4 exp(-a 2 fc 2 /2n 2 ) and S k = (1 - log(47ra 2 ))/2, 
once again displaying a relation between localization in 
momentum space and entropy in that space. It is then 
natural to introduce an entropy in phase space as 



S(i/>) = S x {i>) + S k (i/>) 



(4) 



This definition respects invariance under dilatations (as 
also does AxAp). In the entropic approach the uncer- 
tainty principle is expressed by the fact that S x and 
Sk cannot both be lowered simultaneously without limit. 
More precisely, the entropy in phase space satisfies the 
lower bound jjj El El 



S(V0 > <2(l-log2) 



(5) 



This bound is sharp and is attained by the Gaussian func- 
tions (with non complex width). These are the unique 
mimmizers 0. In the one dimensional case this bound 
implies Heisenberg's relation AxAp > h/2 |fj|, |Ï6|] . We 
note that, in the litcrature, the name phase space entropy 
is used with various meanings, and frequently refers to 
entropies of matrix densities defined on phase space by 
means of a Wigner or a Husimi transformation (Wehrl 
entropy) ||. S can also be viewed as the information 
entropy of a probability density in phase space, namely, 
p( X> k) = \iP(x)\ 2 \ï,(k)\ 2 . 

2. Consider now two partides constrained to have a 
symmetric or an antisymmetric relative wave function 
(say, two electrons in a singlet or triplet spin state, re- 
spectively). One can ask which is the most localized state 
in phase space in both cases. In the symmetric case the 
answer is just a Gaussian function since this is the abso- 
lute minimum of S(ip) and it is symmetric. In the anti- 
symmetric case the solution is not known but certainly 
there exists a sharp lower bound, which will be denoted 
C d 

SW >C d >d(l- log 2) , i;{x) = -V(-aj) ■ (6) 

That the bound Cd is sharp implies that there exists 
a (minimizing) sequence of antisymmetric functions tp n 
such that linin^oo S(ip n ) = Cd- It is not guaranteed, 
however, that S(ip) — Cd is fulfilled for any function 
ip(x) in L 2 (R d ). 

In view of the fact that in the symmetric case the min- 
imizer is a Gaussian, the ground state of a quantum har- 
mònic oscillator, a natural guess would be to take the 
first excited state for the minimizer in the antisymmetric 
case. In one dimension this gives 
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where 7 is Euler's Gamma. Numerically S = 0.847579. 
This is not, however, the true minimizer, since in |fÏ8[ a 
better guess ipo was proposed which yields 



5(^o) = 2(l-log2)>Ci 



(8) 



(numerically S = 0.613706). In fact, based on strong 
numerical evidence, this was conjectured in [pi to be 
the true sharp bound in one dimension. (We note that 
tj^o represents actually a minimizing sequence and not a 
truly acceptable state.) 

3. Whether this conjecture is true or not, one can 
now put an upper bound on Cd for arbitrary d. Let 
ipo(x) be the one-dimensional best guess for the anti- 
symmetric case found in [Q and let g(x) be a (d — 1)- 
dimensional Gaussian function, then for the antisymmet- 
ric d-dimensional function 



i/j(x) = ipo(xi)g(x 2) ■■■ ) Xd) 



(9) 



we can immcdiatcly evaluate the entropy noting that this 
quantity is additive for separable functions 

S(V0 = S(Vo) + S(g) . (10) 

This implies the bounds 

(d + l)(l-log2) >C d >d(l-log2). (11) 

The Fourier transform can be regarded as an opera- 
tor from LP(R d ) into L q (R d ), where p' 1 + qr 1 = 1 and 
1 < P < 2 < q |Ï9| . As already noted in |H| the mini- 
mization of the functional S in R d is directly related to 
the norm of the Fourier transform operator, defined as 
the supremum 

01 HV'llg/llV'llp on Rd - [ This is because 
derivating J d d x\tp\ p or (2nh)~ d f d d k\ip\ q with respect q 
at q = 2 generates S x (t(j) and Sk(tp) respectively.] The 
relation between norm and entropy extends immediately 
for the subspace of antisymmetric functions fÏ8|| . Let 
Kd^q be the norm of the Fourier operator on the anti- 
symmetric subspace of M. d . For this quantity using the 
same separable function introduced in we find the 
bounds 

(pl/p 9 -l/9)(«í+l)/2 < Kd q < (pl/Pg-l/fljl/2 _ (12) 

4. The form of tpo(x) is quite remarkable. It is an an- 
tisymmetric array of very narrow Gaussian functions lo- 
cated at equally spaced points and with amplitudes mod- 
ulatcd by a very wide Gaussian function. Namely, 



-n-\flo? 



, a-> H 
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As indicated, a is a vanishingly small positive parameter. 
(We will implicit assume this in what follows.) Alterna- 
tively one can use 



(14) 



Both expressions coincide for small a. This function is 
normalizcd. In the sense of distributions, ïpo(x)/a is 
equivalent to 



d {x) 



J2(-l) n S(x-n 



(15) 



However do does not have a well defined entropy: the 
localization at any particular x = n + | contributes with 
—00 to the entropy (an ultraviolet divergence) but the 
freedom to choose any n adds a +00 (an infrared diver- 
gence). This refers both to the position and the momen- 
tum space entropies since they are equal for do (do is 
essentially equal to its Fourier transform). 

The function ipo{x) is a regularized version of do in 
the infrared sector (by using e~ 7m x instead of 1) and in 
the ultraviolet sector (by using e~ 7TX l a instead of S(x)). 
Remarkably the actual value of the entropy depends on 
the concrete choice of the regularization, i.e. there is no 
natural value to be assigned to the entropy of do- 

In order to analyze this let us consider a more general 
regularization 



xl>(x) = V (ax)^2(-l) n ^(x-n-í)/a) 

= ^(-l)>(a(n+i)K((x-n-i)/a) 



(16) 



(As always, the equality refers to a — > + .) Hcre ip(x) 
and £(x) are two even smooth functions of ràpid fall at 
infinity which we assume to be normalized. In this case 
if)(x) is also normalized and it will be instructive to show 
this explicitly. Because the £(x) falls rapidly at infinity 
the regularized deltas do not overlap for small a, thus 

ídx\i>(x)\ 2 = /dx^|^(a(n+i))| 2 |C((x-n-I)/a)| 2 

•* «' .„ r n 



^|^(a(n+i))| 2 í dxMx- 
a^l^n+i))! 2 í dx\Ç(x)f 



è)/«)l 2 



dy\tp(y)\ 2 / dx|e(x)| 



= 1. 



(17) 



A similar calculation for the entropy in position space 
gives 



S X W = S x (<p) + S X (Ç) . 



(18) 



In order to compute the momentum space entropy, we 
will use from now units 2?:h = 1 which give simpler for- 
mulas (in particular Sk(ip) = S x (ip)). Note that Sk(ijj) 
has been defined so that it is numerically independent of 
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H (for fixed ip(x)). Using Poisson's summation formula 
J2 n ( t>( n ) = En <K n )> one easüy obtains 



4,{k) = -í|(afc)^(-l)X(fc-n-i)/ a ) 

= -»Xi(-in"(a(n+W-n-5)/a) 



(19) 



This expression is formally idcntical to ip(x) in ( |16[) re- 
placing tp and £ by £ and y> and so 

S k {^) = S X (0) + S x (|) = S fc M + s k (0 . (20) 



This implies for the phase space entropy 
= SfoO + 5(0 . 



(21) 



The minimum is then obtained by choosing Gaussians as 
the regularizing functions tp and £, as in ipo(x). 

5. We will now study maximally localized (minimum 
uncertainty) states in other subspaces of L 2 (R) . To this 
end let us consider generalizations of the distribution 
do(x). These more general distributions u{x) will be com- 
posed of a set of well-separated Dirac deltas with different 
amplitudes and enjoying some periodicity properties. Af- 
ter suitable infrared and ultraviolet regularization we can 
then use the same procedure described above to computo 
the entropy. Specifically, we assume 

N 

u{x) = bkS(x — Xk) for < x < r , 



< x\ < ■ ■ ■ < xm < r , 
\u(x + r)\ = \u(x)\ . 



(22) 



Assuming that no bk vanishes, u(x) is the superposition 
of N series, each series k — 1,...,N being composed of 
equidistant deltas of strength |òfe| 2 . The period r is com- 
mon to the N series. In principle one could consider a 
larger family of distributions by taking linear combina- 
tions of m's, J2í u í( x ) w hh different periods J-j, however, 
in order to guarantee that the deltas are well-separated 
(and this is essential for being able to computo the norm 
and the entropy) the ratios tí/tj must be rational num- 
bers. In this case we are back within the class of distribu- 
tions defined in (|2^) with r equal to a common múltiple 
of the fj. 

Next we regularize u(x) with normalized functions 
<p(x) and £(x) as in (fÏ6|). This amounts to make a con- 
volution of u{x) with Ç(x/a) and multiply by <p{ax) (in 
any order). The norm of the resulting function ij)(x) can 
be computed along the lines of ( |Ï7| ) (working with each 
series of deltas separately). This gives 



/ 



N 



"Ei 



(23) 



The norm decreases for increasing r (for fixed bk) since 
less strength falls under the profile tp(ax). The entropy 
in position space can also be computed (exploiting again 
that the deltas are well-separated as a goes to zero) and 
this gives 

S x (V) = S x (<p) + S x (0 + S(b)- log(r) , (24) 



where 

S ( b ) = -^Pklogpk , Pk 



N 



N 



fe=l 



I6 fe i 2 /Ei^i 2 - ( 25 ) 



Eq. (gj) nicely shows how the various structures com- 
bined in tp( x ) contribute additively to the entropy. S(b) 
corresponds to an entropy for the mixing of the various 
series in u{x). If the number of series is small, or more 
precisely if the strength is concentrated in a few series, 
this mixing entropy will decrease. However, the term 
— log(r) indicates that what really matters for the en- 
tropy of ip(x) is the effective number of series per unit 
length. This is sensible since any u{x) with N series and 
period r can also be regarded as having nN series and 
period nr for n = 2, 3, . . .. 

Assuming that ü(x) belongs to the same class ( ^2f ) for 
some bk, N and r, the total phase space entropy becomes 

S(V0 = S(<p) + S(0 + S(b) + S(b) - log(rr~) . (26) 

As a consequence the optimum small and large scale pro- 
files ^ and tp are Gaussians, and we will assume that in 
what follows. 

We note that similar relations are obtained for the p- 
norms, namely 



dx\ip(x)\ p = r 



N 

-1 \ " 

k=l 



\b k \* / dy\p(y)\P / dx\Ç(x)\> 



6. The Fouricr transform operator naturally decom- 
poses L 2 (R) into the four invariant subspaces correspond- 
ing to its four eigenvalues A = ±l,±i. We want to find 
maximally localized states in each such subspace. The 
Gaussians belong to the subspace A = 1 so they are 
also the minimizer in this subspace. Our best guess for 
antisymmctric functions ipo(x) belongs to the subspace 
A = —i, so this is also our best guess in this subspace. It 
remains to consider the cases A = — 1 and A = i. Brutc 
force numerical minimization suggests that the minimiz- 
ers in these two cases belong to the class (^2|) with Gaus- 
sian regularization. In order to study this further let us 
introduce the set of distributions 



> a ,P) = E' 



S(x — r(n + Oé)) 



(28) 



for non vanishing real r and real a and (3. Not all these 
functions are independent since 

(f>(x,r,a, j3) = 4>{x,—r,—a,—0) = <t>[x,r,a 1 j3 + l) 

= e~ 2mí3 cb(x,r,a + l,(3) (29) 
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and all difïerent results are covered byO<r, < a, (3 < 
1. Undcr Fourier transform onc finds 

4>(x, r, a, (3) = r^e 2 ™^^, r -1 , -0, a) . (30) 

Therefore, the distributions <fi(x, r . a, (3) and their Fourier 
transform belong to the class in ( Bg ) with N = 1. 

In order to minimize the mixing entropy 5(6) we start 
by looking for solutions of the form cj> — X<p. The unique 
solutions are <j>(x, 1, 0, 0) with À = 1 and <j){x, 1, 5, 5) with 
A = — i both with entropy 1 — log 2 after Gaussian regu- 
larization (the latter solution is just do). 

The next simplest guess is to project a single distribu- 
tion <p on the subspaces A = — 1 or À = i. If F denotes 
the Fourier transform operator, the projector on the sub- 
space A is P x = (1 + X^F + \~ 2 F 2 + A- 3 F 3 )/4. Thus 
a combination of no more than four 0's is sufficient to 
yield a distribution of the type ü(x) = Xu(x). A further 
simplification can be achieved as follows. All functions 
with A = — 1 are necessarily even (symmetric), so for 
these functions P\=-i = (1 — F)/2 and only two </>'s are 
involved after projection. Similarly, the functions A = i 
are odd and for these functions P\=i = (1 — iF)/2. In 
both cases we will seek a minimum of this form, i.e. pro- 
jections of a single even or odd cj). Of course it is not 
guaranteed that the true minimum is obtained in this 
way. 

Two remaries should be made. First, the projection 
P\4> involves periods r and r _1 . This combination is ac- 
ceptable only if r/r -1 = r 2 is a rational number. We will 
take r = y/q/p, where the positive integers q and p have 
no common divisor. Second, although each <f> contains a 
single series of deltas, the combination of two <^>'s may 
have many more than two series of deltas. For instance 
u(x) = 4>(x,l/2, 0, 0) + <t>(x,2, 0, 0) has period r = 2 and 
N = 4 with x k = 0,1/2,1,3/2 and b k = 2,1,1,1. In 
practice u(x) will be composed of two <^'s with periods 
r = \J l/p an d r~ 1 . In order to disentangle the different 
series contained in u(x) it will be necessary to bring all 
the <f>'s to a common period. To this end the following 
identity is useful 

n—l 1 
<t>{x, r, a,(3)=J2 e~ 2 ^ k <P{x, nr, ^U/?) , (31) 

fc=0 

where n is any positive integer. In our case the smallest 
common period is r' = pr = qr^ 1 = ^fqp. 

Case A = — 1. There are three classes of single dis- 
tributions 4> which are even under x — > —x, namely 
(/)(x,r,0, 0), 4>(x,r, 0, 1/2) and (/>(x,r, 1/2,0), and arbi- 
trary period r. The case (1/2,0) needs not be consid- 
ered since it will be generated after projection of the case 
(0,1/2). 

Projection of <p(x, r, 0, 0) on the A = — 1 subspace gives 

u(x) = <f>(x,r,0,0) - r^ifrfar- 1 , 0,0) . (32) 

In this case we find a total of N = q+p—1 series of period 
r ' = \/QP- (Actually <p(x, r, 0, 0) gives p series of period 



r' and (p(x, r _1 , 0, 0) gives q series of period r' but one of 
the series is common to both functions.) Of these, p — 1 
series have amplitude bk = 1, q — 1 series have amplitude 
bk = —r^ 1 and there is a single series with amplitude 
1 — r _1 . A numerical survey shows that the best case 
corresponds to q = 1 and p = 2 (or vice versa). This 
corresponds to a distribution u(x) composed of only two 
series of period r' = v2 and weig hts pi = (2- V2)/4 an 
pi = (2 + \/2)/4. The mixing entropy is 5(6) - log(r') = 
log 2 + -^log(V2 — 1). The total phase space entropy 
(i.e., adding the Gaussian contribution 1 — log 2 and an 
overall factor of two to account for the momentum space 
entropy) is thus 

S x =-i = 2 + v^logív^- 1) = 0.753550 (33) 

corresponding to the distribution 

u(x) = (1 - V2)c/)(x, y/2, 0, 0) + ct>(x, y/2, 1/2, 0) . (34) 

This is our best guess for the subspace A = — 1. 

The other possibility in this subspace is to project 
cf)(x, r, 0, 1/2). The analysis is slightly more involved in 
this case. If p is odd the number of series is p + q and 
the total entropy is 5 = 2 for any choice of p and q. If 
p is even the number of series is q + p — 1 and the en- 
tropy depends of the concrete vàlues of q and p. The best 
choice corresponds to q = 1, p — 2 and yields the same 
distribution u{x) in (p4). 

Case A = +i, The only possibility of an odd (j> is 
4>{x,r, 1/2,1/2). This gives 

u{x) = cf>(x, r, 1/2, 1/2) - r-^ix, r" 1 , 1/2, 1/2) . (35) 

If p or q is even the number of series is q + p and for all 
choices of q and p the total entropy is 2. If q and p are 
odd the number of series is q + p — 1 . The best choice 
corresponds to q = 3 and p = 1 (or vice versa) with a 
total phase space entropy 

S x=+i = 2 - -|= log(V3 + 1) = 0.839465 (36) 
v3 

This is our best guess for the subspace X = +i. 

These results are well below the guesses based on the 
lowest harmònic oscillator states which yield 5 = 1.15934 
for A = -1 and 5 = 1.38155 for X — i. 

By using (|27j ) these results extend immediately to rela- 
tions for the norm of the Fourier transform operator from 
L p (R d ) into L q restricted to the proper subspaces 
A = ±l,±i, 
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